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Abstract: We consider stationary solutions of the three dimensional Navier-Stokes 
equations (NS3D) with periodic boundary conditions and driven by an external 
force which might have a deterministic and a random part. The random part of 
the force is white in time and very smooth in space. We investigate smoothness 
properties in space of the stationary solutions. 

Classical technics for studying smoothness of stochastic PDEs do not seem to ap¬ 
ply since global existence of strong solutions is not known. We use the Kolmogorov 
operator and Galerkin approximations. We first assume that the noise has spatial 
regularity of order p in the based Sobolev spaces, in other words that its paths 
are in H^. Then we prove that at each fixed time the law of the stationary solutions 
is supported by 

Then, using a totally different technic, we prove that if the noise has Gevrey 
regularity then at each fixed time, the law of a stationary solution is supported 
by a Gevrey space. Some informations on the Kolmogorov dissipation scale are 
deduced. 

Key words: Stochastic three-dimensional Navier-Stokes equations, invariant mea¬ 
sure, Gevrey spaces, Kolmogorov operator, Kolmogorov dissipation scale. 


Introduction 


We are concerned with the stochastic Navier-Stokes equations in dimension 3 
(NS3D) with periodic boundary conditions and zero mean value. These equations 
describe the time evolution of an incompressible fluid and are given by 


( 0 . 1 ) 


dX + v{-X)X dt + (A, V)A dt + Xpdt = (j)iX)dW + g{X)dt, 

(divA)(t,^) = 0, for ^ G D, t > 0, 

J^X{t,^)d^ = 0, for t > 0, 

A(0,O = xoiO, for^GD, 


where D = (0, 27r)^. We have denoted by X{t,^) the velocity and by p{t,^) the 
pressure at time t and at the point ^ G D, also v denotes the viscosity. The external 
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force acting on the fluid is the sum of a random force of white noise type (j){X)dW 
and a deterministic one g{X)dt. 

As it is well known, in the deterministic case, global existence of weak (in the 
PDE sense) solutions and uniqueness of strong solutions hold for the Navier-Stokes 
equations. In space dimension two, weak solutions are strong and global existence 
and uniqueness follows. Such a result is an open problem in dimension three (see 
m for a survey on these questions). 

In the stochastic case, the situation is similar. However due to the lack of 
uniqueness, we have to work with global weak (in the PDE sense) solutions of the 
martingale problem (see [ 7 ] for a survey on the stochastic case). Roughly speaking, 
this means that in (EU, we take X, p and W for unknown. 

As is usual in the context of the incompressible Navier-Stokes equation, we get 
rid of the pressure thanks to the Leray projector. Let us denote by {X, W) a weak 
(in the PDE sense) stationary solution of the martingale problem (ED and by p 
the the law of X(t ), which is an invariant measure if we can prove that 10.111 defines 
a Markov evolution. In this article, we establish that /i admits a finite moment in 
spaces of smooth functions provided the external force is sufficiently smooth. We 
think that this is an interesting question to study. First, it can be seen that if we 
were able to prove that p has a moment of sufficiently high order in a well chosen 
Sobolev norm (order 4 in or 2 in for instance) then this would imply global 
existence of strong solutions for p, almost every initial data. 

Moreover, this result is an important ingredient if one tries to follow the method 
of (21 to construct a Markov transition semi-group in HP{D) under suitable con¬ 
ditions on (j) and g. Since even uniqueness in law is not known for NS3D, such a 
result might be important. 

We first prove that if the external force is in and the noise term has 

paths in HP{D) then p admits a finite moment in the Sobolev space 

Note that analogous results are well-known for the two dimensional Navier- 
Stokes equations (NS2D). Actually a stronger result is true for NS2D. Namely, for 
any square integrable xq, the unique solution of NS2D is continuous from (0, oo) 
into HP{D) and is square integrable from into It follows that p 

admits moments of any orders in HP{D) and a moment of order 2 in 
This stronger result is linked to the global existence of strong solutions for NS2D. 

This kind of idea cannot be used for NS3D and we use a generalization of an idea 
used in [ 2 ] for the case p = 1. The method is based on the use of the Kolmogorov 
operator applied to suitable Lyapunov functional. These functionals have already 
been used in the deterministic case in dn, chapter 4. 

Using a totally different method, we establish also that the invariant measure 
p admits a moment in a Gevrey class of functions provided the external force has 
the same regularity. Gevrey regularity has been studied in the deterministic case 
in dni and d3- Our method is based on tools developed in dOl- In dH’ these 
tools have been used to obtain an exponential moment for the invariant measure in 
Gevrey norms in the two dimensional case. The arguments used in da, da do not 
generalize to the three dimensional case since there strong existence and uniqueness 
is used. The three dimensional case NS3D requires substantial adaptations. We 
develop a framework which gives a control on a Gevrey norm by using a control of 
the {D)-noim of X only at fixed time. 
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Actually, in this way, we are able to generalize to NS3D the results of d, m 
However, we do not have exponential moments. We deduce that the Kolmogorov 
dissipation scale is larger than . This is certainly not optimal since it is expected 
that the scale is of order 1 / 4 . Note that our result is rigorous and does not use any 
heuristic argument. 


1. Notations 


For TO G N, we denote by the space of functions which are restrictions of 

periodic functions in and whose average is zero on D. We set 

H = {X G I div X = 0 on D } , 

and 

Let TT be the orthogonal projection in L‘^{D)^ onto the space H. We set 
A = tt{-A), D{A) = VnE.l^^{D) and B(m) = TT ((u, V)u) . 

It is convenient to endow with the inner product ((•, ■))m = {A^ •, A^ 

The corresponding norm is denoted by || j|^. It is classical that this defines a norm 
which is equivalent to the usual one. For to = 0 we write | • | = || jig and for to = 1 
we write || j| = || j|j^. Note that, since we work with functions whose average is zero 
on (0, 27r)^, we have the following Poincare type inequality 

W^Wm, ^ II^ILs , rni < m 2 , X G 

We also use the spaces Lp{D)^ endowed with their usual norm denoted by Ijp. 
Moreover, given two Hilbert spaces Ki and K 2 , C 2 {Ki-,K 2 ) is the space of Hilbert- 
Schmidt operators from Ki to K 2 ■ 

The noise is described by a cylindrical Wiener process W defined on a Hilbert 
space U and a mapping (fi defined on H with values in C 2 {U ; H). We also consider a 
deterministic forcing term described by a mapping g from H into H. More precise 
assumptions on cj) and g are made below. 

Now, we can write problem in the form 

{ dX + vAXdt + B{X)dt = (j){X)dW + g{X)dt, 

X(0) = xo. 


In all the paper, we consider a iL-valued stationary solution (X, IT) of the martin¬ 
gale problem jnj. Existence of such a solution has been proved in . We denote 
by /i the law of X(t). We do not consider any stationary solutions but only those 
which are limit in distribution of stationary solutions of Galerkin approximations 
of 11 .111 . More precisely, for any X G N, we denote by Pv the projection of A 
onto the vector space spanned by the first X eigenvalues and consider the following 
approximation of am 


( 1 . 2 ) 


{ dXjsf vAXj^dt P[,fB(Xj^)dt 

x^(0) 


PN(t){XN)dW -I- PNg{XN)dt, 
PnXq. 


It can be easily shown that n~a has a stationary solution Xjy. Proceeding as 
in [^, we can see that their laws are tight in suitable functional spaces, and, up 
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to a subsequence, {Xn,W) converges in law to a stationary solution {X,W) of 
(11.111 . Actually the convergence holds in C(0, t; (H for 

any t, s > 0. We only consider stationary solutions constructed in that way. 

To obtain an estimate for stationary solutions of (II.Ill (limit of Galerkin ap¬ 
proximations), we proceed as follows. We first prove the desired estimate for every 
stationary solutions of and then we take the limit. 

The reason why our results are only stated for solutions limit of Galerkin ap¬ 
proximations comes from the fact that it is not known if computations applied to 
solutions of Galerkin approximations can be applied directly on solutions (A, W) 
of the three-dimensional Navier-Stokes equations. 

Some of our results describe properties of /i in Gevrey type spaces. These spaces 
contain functions with exponentially decaying Fourier coefficients. According to 
the setting given in cni, we set for any (a, /3) G R+ x (0,1] 


I ll^llG(a,/3) 

[ G{a,(3) 


aA 2 , 


As e 




= \xm\ 


\m\G{a,l3) < OO 


} 


where (i(fc))fcez3 are the Fourier coefficients of x G H. Moreover, for any {x,y) G 
G(a,/3)^, we set 


y] = ^ \kf (x{k)m) ■ 


Clearly, \^G{a, ■)G(a 13 ) J ^ Hilbert space. 

We are not interested in large viscosities and in all the article it is assumed that 
1 / < 1. We will use various constants which may depend on some parameter such 
as p, ... When this dependance is important, we make it explicit. 


2. ]HI(!er(D)-regularity 

Let p G N. We now make the following smoothness assumptions on the forcing 
terms. 


Hypothesis 2.1. The mapping (j) (resp. g) takes values in /I 2 (G;iF (H ]HIPg^(Zl)) 
(resp. H(D)) and 4>: H ^ £2 (H;iL n ElPg^(Zl)) and g : H ^ Hr]WPjj}-{D) 
are bounded. 


We set, when Hypothesis o holds, 

(ll^ll£2(C/;Hger(i>)) + llsllp-l) ■ 

It is also convenient to define 

Bp = sup \\(t>\\%^u-,nU(D)) + I sup \\g\\l_^. 

The aim of this section is to establish the following result. 

Theorem 2.2. Let p he the invariant law of a stationary solution X of the three 
dimensional Navier-Stokes equations that is limit of stationary solutions of Galerkin 
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approximations. Assume that Hypothesis lit. 1\ holds for some p > 1. For any < 1, 
there exists Cp^n depending on p, v and Bp such that 

f dfiix) < Cp^^. 

J H 

Let us make few comments. 

Note that it would be very important to obtain an estimate on Ikllp+i dfi{x) 
with pSp > 3. Indeed, by Agmon inequality , we have 

[ \x\lo dn{x) < c [ \x\^-^^P \\x\\i dfi{x) 

JH J H 

and this would give an estimate on the left hand side. Since uniqueness is easily 
shown to hold for solutions in L'^{0,T; L°°{D)^), a classical argument could be 
used to deduce that for fi almost every initial data there exists a unique global 
weak solution. Combining with the result in |^, this would partially solve Leray’s 
conjecture. 

Consider the case g = 0, U = H and (p = Then Hypothesis 12 . II holds for 

any p < s and the unique invariant measure of the three dimensional linear stochas¬ 
tic Stokes equations in H is in ]HIp+^(D) with probability zero if r > s. Therefore 
it seems that IHIp.i.i is the strongest norm we can control under Hvpothesis l2.1l 

Remark that in the two dimensional case a much stronger result holds. Indeed, 
standard arguments imply that under Hypothesis 12.11 we have for any invariant 
measure fi and any q G N* 

[ \\x\\p‘ d^i{x) < oo, f ||x||p+i d/r(x) < oo. 

J H J H 

In the proof, we use ideas developped in CHI Similar but more refined techniques 
have been used in CD to derive interesting properties on the decay of the Fourier 
spectrum of smooth solutions of the deterministic Navier-Stokes equations. Using 
such techniques does not seem to yield great improvement of our result. Indeed, 
trying to do so, we have been able to improve the estimate of Theorem l2.2l as follows 


[ ll^llpii d^x) < 2Bp + c2P(l + Bo), 

JH 

where c and c* are positive constants and c* is close to 1.02. We have not been able 
to derive very interesting results from this improved estimate and therefore have 
preferred to give the simpler one which follows from easier arguments. 

Proof: The proof is rather standard. We do not give the details. 

Let {nN)Nen be a sequence of invariant measures of stationary solutions (X^)]^ 
of (11.211 such that there exists a subsequence {Nk)kGfi such that converges to X 
in law. It follows that (/i 7 Vj,)fegN converges to p, (considered as probability measures 
on D{A~^)). 

We denote by Ljv the Kolmogorov operator associated to the Galerkin approxima¬ 
tion of the stochastic Navier-Stokes equations 

LNfix) = ^tr {{PN<p){x){PN(p)*ix)D^f{x)) - {uAx + B{x) - g{x), Df{x )), 


for any / € C'^{PmH-,M.) and x G PnH. 
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The proof of Theorem is based on the fact that, for any N £N, we have 


(2.1) 


/ LNf{x)dgN{x) = 0, 




J PnH 


provided f € 

(PnH 

;R) verifies 



*) 

IPnH I/(^)I dgN^x) 

< oo 

(2.2) 

ii) 

Jp^H \LNf{x)\dgN{x) 

< oo 


in) 

.. 

Jp^H \ iPN(j))*{x)Df{x)\^ dflN{x) 

< oo 


It follows from [ 7 ], Chapter 1.2, Corollary 1.12 that for any p S N 
(2.3) J =e(|Xjv(0)|"p) < 2 ^^,, < oo. 


and 

(2.4) 


'PnH 


' d^N{x) < Bq. 


The result in [ 7 ] is given for g = 0 but the generalization is easy. 
Thanks to , we use EH with 


/ = 




£r) — 


2p-l 


We obtain after lengthy but easy computations and some estimates on the nonlinear 
term borrowed from CHI, chapter 4 (in particular equation (4.8)) that there exists 


Cp such that 

(2.5) 

with 


Bp ‘^Bp + CpBQ + 1 , 


Rp = y 


1+lla^l 


p+i 


Then arguing as in ca, chapter 4, we set 


l+£'i 


-dnN{x). 


Mp = V 


f , ,n1 /2 p-1 

/ (l+l|a^|U dgi^{x), 

J PnH ^ ' 


and deduce 

Mppi < 

which yields 


(2.6) 

[ 

JPnH 

It is then standard to deduce the re 


a;||p+Y d[ij^{x) < Cp^^. 


Fatou Lemma and lower semicontinuity of 
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3. Gevrey regularity 

3.1. Statement of the result. We now state and prove the main result of this 
paper. It states that if the external force is bounded in a Gevrey class of functions, 
then ^ have support in another Gevrey class of function. 

The main assumption in this section is the following 

Hypothesis 3.1. There exists {a,P) £ x (0,1] such that the mappings g : H ^ 
G{a,P) and 4> ■. H ^ C 2 {U;G{a, (3)) are bounded. 

We set 

B'o = sup U{x)\\l,(u,G(a,fi)) + sup ll 5 (a:)llG(a./ 3 ) ■ 

x£H x£H 

The aim of this section is to establish the following results proved in the following 
subsections. 

Theorem 3.2. Let p, be the invariant law of a stationary solution X of the three di¬ 
mensional Navier-Stokes equations that is limit of stationary solutions of Galerkin 
approximations. Assume that Hypothesis I.V.il holds. There exist a family of eon- 
stants only depending on {a,l3,BQ) and a family {av)u^{p.i) of measur¬ 

able mappings H —> (0, a) such that for any v £ (0,1) 

(3.1) f + 

J H 

(3.2) / {a^{x))~^ dfi{x) < Kj{l -\- Bo)i^~X 

J H 

for any 7 £ (0,1). 

This result gives some informations on the Kolmogorov dissipation scale. Indeed, 
it follows from ijOl . 1221 ) that 

where {x{k))k^i3 are the Fourier coefficients of x. 

Hence, if Hypothesis 10 holds with (3=1 and 5 = 0, then |a;(fc)| decreases faster 
than any powers of |fc| for |A:| » (va,^{x))~^. By (13.211 . for any d > 0 

—\-T < C 5 with / cs j/(x) ^u+u p,(dx) < 85 < 00 , 

a^[x) ’ J 

and 05 not depending on u. It follows that |ai(fc)| decreases faster than any powers 
of |fc| for |fc| » ^“(6+5'*). This indicates that the 3D-Kolmogorov dissipation 
scale is larger than Note that by physical arguments it is expected that the 

3D-Kolmogorov dissipation scale is of order of uJ. 

Making analogous computations, we obtain that, for 5 7 ^ 0 and (3 £ (0,1], the 
3D-Kolmogorov dissipation scale is larger than for d > 0. 

In ^1], ^7], analogous results to Theorem 13.21 are proved for the 2D-Navier- 
Stokes equations. Moreover, it is deduced in m that the 2D-Kolmogorov dissi¬ 
pation scale is larger than . In it is established that the 2D-Kolmogorov 
dissipation scale is larger than 1 / 2 , which is the physically expected result. 

We can also control a moment of a fixed Gevrey norm. 
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Corollary 3.3. Under the same assumptions, there exists a family {Cj^a>,i 3 ',u)'y,a>,f 3 ',u 
only depending on {a, P, Bq,v) such that 

(3-3) J (ln+ dp{x) < 

where In'*" r = max{0,lnr} and provided a' > 0 and /3',7 > 0 verify 

P' < P and 2j < — 1. 

3.2. Estimate of blow-up time in Gevrey spaces. 

Before proving Theorem 13.21 we establish the following result which implies that 
the time of blow-up of the solution in Gevrey spaces admits a negative moment. 


Lemma 3.4. Assume that Hunothesis \S. 1\ holds. For any N, any stationary so¬ 
lution Xn of and any v G (0,1), there exist K only depending on (a,/3, ilg) 

and a stopping time > 0 such that 


E sup ||^iv(t)||G(i.t,/3) 
Vie( 0 ,r«) 


< 


-{Bq 1), 


(3.4) 

(3.5) 

This result is a refinement of the result developed by Foias and Temam in cni and 
is strongly based on the tools developed in this latter article. We denote hy p,N the 
invariant law associated to X^. Let us set 


'{T^<t) < K{Bo + l)tK- 


(3.6) 

Clearly 


T^ = 


inf I 


t > 0 


1 -f 


ll^^(i)llG(.q/3) >4(ll^iv(0)f-hi)}. 


E 


^^SUP^^ ll^^(t)llG(.q/5)] < 4E (||X^(0)f + 1 


and ea follows from (EH). It remains to prove (13.511 . 
We apply Ito formula to |l-^Ar(i)|lG(,,t for t G (0,a) 

2 


(3.7) 


d||^iv(i)llG(.q/3) + 2v\\A2XN{t) 


G{vt,l3) 


dt = 


A^iXNp) 


dt 


G{vt,l5) 

-\-dAI(t^ Ilt'jdt, 

where 

m = 2Ig{t) + 2lB{t)+U{t), lB{t) = -(X^(t),i3(X^(t)))c(,*,^), 

l4,{t) = \\PN4>{XN{t))\\'^B 2(U-,G(vt,l3)) ’ h (t) = (5(^Af(^)))-^Af(^))G(!.t,/3) ’ 


M{t) = 2 {XN{s),PiXN{s))dW{s))a^,.t,p)- 

do 

The following inequality is proved in |1 Oj for P < 1 

2 c 

H—^ 

G{vt,l3) 




11^^(011 


(3.8) 2/fj(t) < V A^X]\j{t) 

By Hypothesis 13.II we have 

(3.9) I^{t) 2Ig{t) < Bq-\-1. 


G{ut,l3) ■ 
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Combining ll^ . (CT and (E3, we obtain since (3,v <1 

(3.10) d||^iv(t)|lG(.q/3) < dMit) + ^ ||X^(t)fG(.q/3) dt + {Bq + l)dt. 

Applying Ito formula to ^1 + || Arjv(t) , we then deduce from (13.1011 and 

from Hypothesis EH that for any t G (0, a) and any v <\ 


(3.11) 


-d(l+||A^(t)||^(,,_^))'' < dM[t) + Cqu ^ 


dt, 


where Cq = c{Bq + 1) and 


Mit) = 4 r (l + ||A^(s)||G(.t,;5)) ' (X^is), <)>(A^(s))dlD(s))G(.t.0) ■ 

J 0 


Setting 


Tq = inf t S (0, a) 


A 


M{t) > 


4(l+||X„(0)f) 


T ( , 


" \^4Co(l + ||X,v(0 )f)"y ’ 
we obtain by integration of 113.111! on [0, t] for t G (0, t^) 

l + ll^iv(i)llG(.q/3)<4(l + l!^iv(0)f). 
We deduce that and 

(3.12) P < t) < P < i) + P + ||A^(0)|f 

Since /r is the law of Xn{0), we have 


> 


4:C(]t 


p ((i + l|v„(0)f)’ > jA) = e ii, 1 + ll.f > . 

Applying Chebyshev inequality, we deduce from (EH) 

p((l + l|X„(0)llf >j^) 

(3.13) < (Cot)^ / {l+\\xf)d^iN{x) 


< 2iy-i (l + Bo) {Cot)K 

Moreover 

P(ro"<t) =p(^4(^(l + ||A^(0)f)%up,g[o.tAx-]-M(s)) > l) 
< 4E |^(l + ||AA,(0)f) sup^g[0.iAro"] X(s)^ • 


Taking conditional expectation with respect to the cr-algebra J-q generated by 
AAr(O) inside the expectation, it follows 


'(ro^<t) <4E (1 


(c 


||A^(0)f)'E 
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By Burkholder-Davis-Gundy inequality (see Theorem 3.28 page 166 in m) we 
obtain 


E I sup M{s) 

y s£[0,tAT((''] 




0 <cE((M)MtAro'^)| J^o) , 


and 


(ro^<t) <4E[(l + ||X^(0)f)'E((M)MtAro^)| .Fo)]. 


(3.14) 

We have 

(M)(t)=4 /* (l+||X^(s)|lG(.q0)) 


HXn{s)) X^(s) 


ds. 


Therefore 


(A4) (t) < 4 / ^1 + 
0 


It follows that 


(M)(tAro^)< 


43 (l + ||X^(0)|f 


Hence we infer from and from H-H < 


llG(i/t./3) 


that 


(3.15) P (to"^ < t) < V^. 

Combining and we deduce 


3.3. Proof of Theorem 13.2L 

Let (/iAr)AreN be a sequence of invariant measures of stationary solutions (Xtv)^ of 
C3 such that there exists a subsequence {Nk)k&n such that Xn^ converges to X 
in law. It follows that (/iArj,)fegN converges to /r (considered as probability measures 
on D{A~^)). 


Setting 


a,y{x) = inf 


|s > 0 


ll^ll 


2 

G(!vs,/3) 



it follows that for any 7 G (0,1) 


(3.16) 



g{,uo.AAA) - 



{ay{x)) = dn{x). 


Hence inj is consequence of (lint . Then in order to establish Theorem lim it is 
sufficient to prove (lint . 

Clearly 


P 




2 

G{ut,i3) 


> 




(Ho + 1)) < P (sup,g[o,^W] ||WAr(s)||e(^ 


*./3) 


> 


4 


(Ho + l)) 


+P (r^ < t) 


where has been defined in section 3.2. Applying Chebyshev inequality, we infer 
from Lemma EH and from the fact that, for any t > 0, /r^r is the law of Xjsfit) 

(3.17) ^iN{Ot)=¥{XN{t) G Ot) < [K + l)il + Bo)tii^-X 
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where 

Ot = |a; G D{A {Bq + l)| • 

Notice that O is an open subset of D{A~^). Hence, since fiN,, M (considered as 
probability measures on D{A~^)), then we deduce from (13.1711 that 

(3.18) < liminf < {K + 1)(1 + 

k 

Notice that 

{cc G D{A ~^), au{x) < t} C Ot, 
which yields, by (imii and /r(i7) = 1, 

(3.19) ^j.{x G H , ai,(x) <t)<{K+ 1)(1 + Bo)t^v~^ ■ 

It is well-known that (ITBll for any t > 0 implies (IQ) . which yields Theorem E3 


3.4. Proof of Corollary 13.3L 

To deduce Corollary 13.31 from Theorem it is sufficient to prove that for any 
{a', a, P', P) G (0, oo)^ x (0,1]^ such that P' < /3, we have 

(3.20) lkllG(a',/3') < exp {c{P,P') {a')^ lkllG(a,/3) • 

Indeed, implies that for any 7 G K.^ 


+11 II 2 

^\\G{yoLu{x),^) 


l|a;|lG(a',/3'))^ < C7 (^c(/3,/3') -k (a')^ {va^{x)) + (in 

which yields Corollary 13.31 provided Theorem 13.21 is true. 

We now establish (IT^ . It follows from arithmetic-geometric inequality that for 
any k G iP 

(3.21) a'|A:|^ < c(/3, P') (a') (a) p-p' -|-a|fc|^. 

Recalling that 


lG(a',/3') 


X! 1^1^ exp {2a' \kf ) |:E(fc)|^ , 


fcez3 


we infer (13.2011 from (TOTli . 
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